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ABSTRACT 

The combination of the first-year Wilkinson Microwave Anisotropy Probe 
(WMAP) data with other finer scale cosmic microwave background (CMB) ex- 
periments (CBI and ACBAR) and two structure formation measurements (2dF- 
GRS and Lyman a forest) suggest a ACDM cosmological model with a running 
spectral power index of primordial density fluctuations. Motivated by this new 
result on the index of primordial power spectrum, we present the first study on 
the predicted lensing probabilities of image separation in a spatially flat ACDM 
model with a running spectral index (RSI- ACDM model). It is shown that the 
RSI-ACDM model suppress the predicted lensing probabilities on small splitting 
angles of less than about 4" compared with that of standard power-law ACDM 
(PL-ACDM) model. 

Subject headings: cosmology:theory — dark matter — galaxies:halos — gravitational 
lensing — early universe — large-scale structure 
PACS numbers:98.80.Es,98.62.Sb,95.35.+d,98.80.Cq 

1. Introduction 



Mapping the mass distribution of matter in the universe has been a major challeng 
modern observational cosmology. One of the direct procedures to weigh matter in th 
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universe is measuring its deflection of light by gravity. The statistics of gravitational lensing 
can provide us with a very powerful probe of the mass distribution of the Universe. It is well 
known that the Jodrell-Bank VLA Astrometric Survey (JVAS) and the Cosmic Lens All- 
Sky Survey (CLASS) (Browne & Myers 2000; Myers et al. 2003; Browne et al. 2003) have 
provided us with observations of strong lensing probabilities for small image separations 
ranging from 0.3" to 3". By comparing predicted lensing probabilities with observations, we 
can examine the mass distributions of dark matter halos, in particular, their inner density 
slopes (Huterer & Ma 2004; Ricotti 2003). 

Based on the Cold Dark Matter (CDM) model, which has become the standard theory of 
cosmic structure formation, the lensing probabilities strongly depend on the density profiles 
of CDM halos. The lensing model is usually described by a singular isothermal sphere (SIS), 
the Navarro- Frenk- White (NFW) model (Navarro et al. 1996, 1997), or generalized NFW 
(GNFW) density profiles of dark halos (Zhao 1996). Li & Ostriker (2002) employed a semi- 
analytical approach to analyze the gravitational lensing of remote quasars by foreground dark 
halos and checked the plausibility of various lensing models. They found that no model can 
completely explain the current observations: the SIS models predict too many lenses with 
large splitting angles, while the NFW models predict too few small splitting angles. They 
therefore further developed a two-population halo model for lensing: small mass halos with a 
steep inner density slope and large mass halos with a shallow inner density slope, concluding 
that a combination of SIS and NFW halo models can reproduce the current observations 
reasonably well. Unlike previous work that directly models the density profiles of dark 
matter halos semi-analytically, Zhang (2004) generalized the density profiles of dark matter 
halos from high-resolution N-body simulations by means of generalized Navarro-Frenk- White 
(GNFW) models of three populations with slopes, a, of about -1.5, -1.3 and -1.1 for galaxies, 
groups and clusters, respectively. He presented the calculations of lensing probabilities using 
these GNFW profiles for three populations in various spatially flat cosmological models with 
a cosmological constant A. He showed that the compound model of density profiles does not 
match well with the lensing probabilities derived from the combined data of JVAS/CLASS. 
Recently, Metcalf (2004) compared predictions on small scale structure for the ACDM model 
by numerical simulations with observed flux ratios, and found that the disagreements between 
monochromatic flux ratios and simple lens models can be explained without any substructure 
in the dark matter halos of primary lenses. However, spectroscopic lensing observations of 
Q2237+0305 require more small mass dark halos than that expected in the ACDM model. 

In particular, a power spectrum of primordial fluctuation, P p (k), should be assumed 
in advance in the calculation of lensing probabilities. Inflationary models predict a ap- 
proximately scale-invariant power spectra for primordial density (scalar metric) fluctuation, 
P p {k) oc k n with index n = 1 (Guth & Pi 1982; Bardeen et al. 1983). The combination of 
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the first-year Wilkinson Microwave Anisotropy Probe (WMAP) data with other finer scale 
cosmic background (CMB) experiments (Cosmic Background Imager [CBI], Arcminute Cos- 
mology Bolometer Array Receiver [ACBAR]) and two observations of large-scale structure 
(the Anglo-Australian Telescope Two-Degree Field Galaxy Redshift Survey [2dFGRS] and 
Lyman a forest) favor a ACDM cosmological model with a running index of the primordial 
power spectrum (RSI- ACDM), while the WMAP data alone still suggest a best-fit standard 
power-law ACDM model with the spectral index ofn«l (PL-ACDM) (Spergel et al. 2003; 
Peiris et al. 2003). However, there still exist the intriguing discrepancies between theoretical 
predictions and observations on both the largest and smallest scales. While the emergence 
of a running spectral index may improve problems on small scales, there remain a possible 
discrepancy on the largest angular scales. It is particularly noted that the running spec- 
tral index model suppress significantly the power amplitude of fluctuations on small scales 
(Spergel et al. 2003; Yoshida et al. 2003). This imply a reduction of the amount of sub- 
structure within galactic halos (Zentner & Bullock 2002). Yoshida et al. (2003) studied 
early structure formation in a RSI-ACDM universe using high-resolution cosmological N- 
body/hydrodynamic simulations. They showed that the reduced small-scale power in the 
RSI-ACDM model causes a considerable delay in the formation epoch of low-mass minihalos 
(~ 1O 6 M ) compared with the PL-ACDM model, although early structure still forms hierar- 
chically in the RSI-ACDM model. Thus the running index probably affect the strong lensing 
process of distant sources by intervening dark matter halos because the lensing probabilities 
strongly depend on the abundance of dark halos formed in the evolution of the universe. In 
this letter, we will present the first calculation of lensing probabilities in a RSI-ACDM model 
to explore the effect of running spectral index of primordial fluctuation on strong lensing 
probabilities. We adopt a spatially flat ACDM cosmological model through this letter. 

The reminder of this paper is organized as follows. We describe strong lensing proba- 
bilities and mass function of dark halos in Section 2 and 3 respectively. The cross-sections 
for producing multiple images are described in Section 4. Our results is shown in Section 5, 
while a conclusion and discussion are given in Section 6. 

2. Strong Lensing Probability 

The lensing probability with image separation larger than A9 is given by Schneider et al. 
(1992) 

P(> Ad) = J p(z s )dz s £ dD J^ dz £ n{M, z)a(M, z)dM. (1) 
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In this expression, D p (z) = c/H /* ^/(1+2)E(^) is the proper distance from the observer to 
the lens at redshift z where the expansion rate of the universe E(z) = ^/f2 m (l + z) 3 + £7a in 
a spatially flat cosmological model, and p(z s ) is the redshift distribution of distant sources. 
The physical number density n(M, z) of virialized dark halos of masses between M and 
M + dM is expressed as n(M, z) = n(M, z)(l + z) 3 and cr(M, z) is the cross-section defined 
in the lens plane for forming multiple images. 



3. Mass Function of Dark Halos 

In the standard hierarchical theory of structure formation, the comoving number density 
of virialized dark halos per unit mass M at redshift z can be given by the Press and Schechter 
(PS) formula (Press & Schechter 1974): n(M,z) = dN/dM = p f(M,z)/M where p is the 
mean mass density of the universe today and, instead of PS formula in this letter, the mass 
function f(M, z) takes the form of an empirical fit from high-resolution simulation (Jenkins 
et al. 2001) 

, 0.301 din A _1 (M,z) , , . . ,, S8 , 

f(M,z) = — _A_^exp(-|lnA- 1 (M,^) + 0.64| 3 - 88 ). (2) 

Here, A(M, z) = A(M)D(z) and D(z) = e(Q(z))/e(Q m )(l + z) is the linear growth function 
of density perturbation (Carroll et al. 1992), in which e(x) = 2.5x/(l/70 + 209x/140 - 
x 2 /140 + x 4 / 7 ) and Q(z) = O m (l + z) 3 /E 2 (z). The present variance of the fluctuations within 
a sphere containing a mass M can be expressed as A 2 (M) = ^ J °° P(k)W 2 (krM)k 2 dk ) 
where W(kru) = 3[sin(A;rM)/(A;^M) 3 — cos(kryi) / '(^m) 2 ] is the Top-hat window function in 
Fourier space and tm = (3M/47rp ) 1 ^ 3 - The power spectrum of CDM density fluctuations 
is P(k) = P p (k)T 2 (k) where the matter transfer function T(k) is given by Eisenstein & 
Hu (1999), and P p (k) is the primordial power spectrum of density fluctuation. The scale- 
invariant primordial power spectrum in the PL-ACDM model is given by P p (k) = Ak na with 
index n s =l and that in the RSI-ACDM model is assumed to be P p (k) = P(ko)(k/k ) n3( - k \ 
where the index n s (k) is a function of length scale 

The pivot scale A; =0.05 h Mpc -1 , n s (/c )=0.93, and dn s /d\n /c=-0.03 are the best-fit values 
to the combination data of the recent CMB experiments and two other large-scale structure 
observations (Spergel et al. 2003). For both PL-ACDM and RSI-ACDM models, the ampli- 
tude of primordial power spectrum, A and P(k ), are normalized to as = A(rM = 8/i _1 Mpc), 
which is the rms mass fluctuations when present universe is smoothed using a window func- 
tion on a scale of 8/i _1 Mpc. 
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4. The Cross-sections for Producing Multiple Images 

The cross-section for producing multiple images relies on the density profile of dark 
matter halos. Based on previous work (Li & Ostriker 2002; Zhang 2004), we model dark halos 
as an SIS for M < M cl and NFW profile for M > M cl , respectively, where M cl ~ 1O 13 /i _1 M 
corresponds to the cooling mass scale (Porciani & Madau 2000; Kochanek & White 2001). 



4.1. The Cross-section for SISs 

The mass density for a lens of SIS is p(r) = a1/2itGr 2 where a v is the velocity dispersion 
(Schneider et al. 1992). Thus the surface mass density of the SIS satisfies £(£) = al/2G^ 
where £ = £ is the position vector in the lens plane. By defining the length scales in the 
lens plane and the source plane as £ — ^{^v/cfdid^g/dg , r/ = ^odg/d^ respectively, we 
can simplify the lensing equation and obtain the image separation of SIS lensing 

A6 = 1.27'(^)Mjf £(*) 2 / 3 , (4) 

where M15 = Mj (10 15 h~ 1 M Q ) is dimensionless halo mass. Here dg and d^ are the angular 
diameter distances from the observer to the source and to the lens object respectively, while 
d^g is the same quantity but from the lens to the source object. Therefore, the cross-section 
defined in the lens plane for forming two images by an SIS lens with splitting angle A9 > A9 
can be expressed as <r(M, z) = it£q$(A9 — A9 ) where d is a step function. 



4.2. The Cross-section for NFW Models 

The GNFW density profile can be expressed in the form p(r) = p s r^/r a (r + r s ) 3 ~ a (Zhao 
1996) where < a < 3. The GNFW density profile reduces to the case of NFW if a — 1. 
The mass of a dark halo within r 2 oo can be defined as M = 4ir J^ 200 pr 2 dr = 47rp s rf /(ci), and 
r 2 oo is the radius of a sphere around a dark halo within which the average mass density is 200 
times the critical mass density of the universe. The function f{c 1 )= f Q Cl x 2 dx/x a {l + xf~ a 
and C\ = r 2 oo/ r s is the concentration parameter and takes the form (Oguri et al. 2001, 2002; 
Oguri 2004) 

Cl (M 15 , z) = c norm ^^[10M 15 ]-°- 13 , (5) 
1 + z 

where z is the redshift of halo, c norm = 8 (Bullock et al. 2001). So p s and r s can be related 
to mass M 15 and redshift z by 
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200 c? 1.626 m! /3 



p. = !{ci( ^ z) y r, = — ^ fc-Mpc, (6) 

where p Co is the critical mass density of the universe today. The lensing equation for the 
GNFW profile is given by y = x — fj, s g(x)/x; £ = xr s and if = yrgd^/d^ are the position vec- 



tors in the lens plane and the source plane respectively, g(x) = J* udu J °° {u 2 + z 2 ) a {u 2 + z 2 ) + 1 
and 

(is = ^f 1 = ^^)(TJ^)(^-h (7) 

L cr p co Ih x Mpc c/Hq 

where fj, s is a parameter on which the efficiency of producing multiple images is strongly 
dependent. Here E cr = is the critical surface mass density and = dfd^ s /dg. 

The lensing equation curves are symmetrical with respect to the origin. Multiple images 
can be formed when \y\ < y cr , where y CT is the maximum value of y when x < or the 
minimum value for x > 0. Generally speaking, there exist three images for \y\ < y CT . 
We will just consider the outermost two images stretched by the splitting angle A# when 
more than two images are formed. Therefore, we can write the cross-section as o (M, z) 
iry 2 r r 2 $ (A9 — A9 ) with A# > A# in the lens plane for multiple images produced by a 
GNFW lens at z. Here splitting angle A6 is given by A6 = r s Ax/df xs 2x r s /df and x is 
the positive root of the lensing equation y(x) = 0. 
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5. Numerical Results 

In this letter, we assume spatially flat ACDM models characterized by the matter density 
parameter Q m , vacuum energy density parameter Qa- For both PL- ACDM and RSI- ACDM 
models, we take cosmological parameters to be the new result from the WMAP: Hubble 
constant h = 0.71, Q m = 0.27, a 8 = 0.84 (Bennett et al. 2003; Spergel et al. 2003) and 
M cl = 1.5 x 1O 13 /i -1 M . As mentioned above, the lensing probability depend strongly on 
the abundance of viralized dark halos. The mass function of dark halos directly involve the 
calculation of primordial power of density fluctuation. According to Eq.(l), it is clear that 
the effect of running spectral index on mass function of halos cause the difference of lensing 
probabilities between the two models. 

Although the redshift distribution of quasars in the JVAS/CLASS survey is still poorly 
known, the prediction of Dunlop & Peacock (1990) model and the CLASS lensing sub-sample 
redshift measurements suggest that the redshift distribution for CLASS unlensed sources can 
be modelled by a Gaussian distribution with mean redshift < z s >=1.27 (Marlow et al. 2000) 
and dispersion a 2 =0.95 (Chae et al. 2002; Browne et al. 2003; Myers et al. 2003). Thus in 
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this letter we adopt the Gaussian redshift distribution of quasars with the mean redshift 
< z s >=1.27 and dispersion cr 2 =0.95. In the definition of cross-section for forming two 
images, we just consider the criterion A8, and neglect another one q r that is the flux density 
(brightness) ratio of multiple images (Schneider et al. 1992). In order to investigate the effect 
of central black holes or bulges on lensing probability, Chen (2003a,b) introduced q r into the 
calculation for lensing cross-section. Due to the existence of central black holes or galactic 
bulges, y cr becomes extremely large when \x\ approaches zero. Thus y cr can be determined 
by the consideration of q r together with A9. However for GNFW halo models in the absence 
of central black holes or galactic bulges, the lensing equation curves are so smooth that we 
do not need to define cross-section by q r . As for SISs, it is necessary to introduce q r because 
both the predicted lensing probabilities and the determination of the cosmic equation of state 
u are quite sensitive to q r (Chen 2004a,b). Our objective in this letter is to examine the effect 
of running spectral index on the predicted lensing probabilities, so for explicit we neglect 
this selection criterion q r for forming multiple images. In such cases, we calculate lensing 
probabilities with image separations greater than A8 according to the compound model 
of SIS and NFW halo profile in both PL-ACDM and RSI-ACDM models. Our numerical 
results are shown in Fig.l together with the observational one from JVAS/CLASS. As we 
expect, there is slight difference between the two models on large image separations, while 
this difference enlarges with the decrease of the splitting angle. More specifically, the RSI- 
ACDM model can reduce the predicted lensing probabilities on small splitting angle of less 
than about 4" compared with that of PL-ACDM model. 



6. Conclusions and Discussion 



Motivated by the new result on the index of primordial power spectrum from a combina- 
tion of WMAP data with other finer scale CMB experiments and other large-scale structure 
observations, we present the first study on the predicted lensing probability in a ACDM 
model with a running spectral index. In a popular cosmological model from the new fit 
results mentioned above: h = 0.71, Q m = 0.27, a 8 = 0.84, we calculate lensing probabilities 
with image separations greater than A8 according to the compound model of SIS and NFW 
halo profile in both PL-ACDM and RSI-ACDM models. From the analysis above, we can 
see that the running spectral index mainly affect the predicted lensing probability on small 
image separation. It is well known that structures in the universe forms hierarchically in 
standard CDM models. Yoshida et al. (2003) found that although this hierarchical forma- 
tion mechanism do not work well in RSI-ACDM model compared with that in PL-ACDM 
model and it also is not clear that the PS theory can be used in RSI-ACDM model, the 
mass function measured by high-resolution cosmological N-body/hydrodynamic simulations 
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overall match the PS mass function for both RSI-ACDM and PL-ACDM model. 

In addition, because the running spectral index model predicts a significant lower power 
of density fluctuation on small scales than the standard PL-ACDM model (Spergel et al. 
2003; Yoshida et al. 2003), it should also attract considerable attention in studies on weak 
lensing by large-scale structure (Ishak et al. 2003), especially on skewness (Pen et al. 2003; 
Zhang et al. 2003) which characterizes the non-Gaussian property of k, field in the nonlinear 
regime. 
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Fig. 1. — The lensing probabilities for image separations greater than A9. The solid-line 
histogram represents observed lensing probabilities from JVAS/CLASS. The solid line is the 
lensing probabilities for running spectral index ACDM model, while the dashed line is that 
for power law ACDM model. 
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